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ABSTRACT 

We report subarcsecond-resolution VLA and Keck mid-infrared imaging of the 
dwarf starburst galaxy II Zw 40. II Zw 40 contains a bright, compact thermal 
radio and infrared source with all the the characteristics of a collection of dense 
HII regions ionized by at least 14,000 O stars. The supernebula is revealed to 
be multiple sources within an envelope of weaker emission. The radio emission 
is dominated by free-free emission at 2 cm, and the spectrum of this emission 
appears to be rising. This suggests that the free-free emission is optically thick 
at 2 cm, and that the individual HII regions are ~1 pc in size. This complex 
of "supernebulae" dominates the total infrared luminosity of II Zw 40, although 
the radio source is less than ~150 pc in diameter. Multiple super star clusters 
appear to be forming here, the much larger analogues of large Galactic HII region 
complexes. 

Subject headings: galaxies: individual (II Zw 40) — galaxies: starburst 

- galaxies: star clusters — galaxies: dwarf — galaxies: peculiar — radio 
continuum: galaxies-infrared 

1. Introduction: Infrared and Radio Supernebulae in Starbursts 

Starburst galaxies are distinguished by the concentration, and not simply the quantity, 
of star formation. Optical and UV studies with the resolution of the HST have found 
star formation regions to contain super-star clusters (SSCs), as populous and as dense as 
globular clusters but with ages of only a few million years. Observations in the radio and 
infrared are finding the even younger and obscured counterparts of these star clusters. 
Traced by optically thick free-free emission at centimeter wavelengths, radio/infrared 
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"supernebulae" (Turner, Ho, & Beck 1998; Turner, Beck, & Ho 2000; Gorjian, Turner, 
& Beck 2001) or "ultra-dense HII regions" (Kobulnicky & Johnson 1999) have been 
discovered in nearby starburst galaxies. The discovery of these regions has been due to 
the development of subarcsecond imaging in the infrared and radio. These HII regions 
are bright and compact, excited by immense clusters of hot young stars, or "starforming 
clumps" (Beck, Turner, & Gorjian 2001). The sources have in common all of the signs 
of young star formation: high gas density (often > 10 4 cm" 3 : Turner, Ho, & Beck 1998; 
Kobulnicky & Johnson 1999; Tarchi et al. 2000; Meier, Turner, & Beck 2002), a cm- wave 
continuum with spectral index from -0.1 to +1.5, substantial extinctions (Av >~ 10 mag: 
Kawara, Nishida, & Phillips 1989; Ho, Beck, & Turner 1990), and very strong mid-infrared 
emission (Gorjian, Turner, & Beck 2001; Beck, Turner, & Gorjian 2001; Dale et al. 2001; 
Vacca, Johnson, & Conti 2002). It is not just that these supernebulae are very luminous in 
the infrared, they may actually dominate the infrared emission of the entire galaxy. In the 
nearby dwarf galaxy NGC 5253, for example, at least 75% of the mid-infrared emission seen 
by IRAS and 25% of the total bolometric luminosity appear to come from one obscured star 
cluster less than 2 pc in diameter (Gorjian, Turner, & Beck 2001). In another starbursting 
dwarf, He 2-10, essentially all the mid-infrared emission is from star-forming clumps in a 
200 pc by 50 pc disk (Beck, Turner, & Gorjian 2001; Vacca, Johnson, & Conti 2002). If 
infrared emission traces and quantifies star formation, as has been thought since IRAS, then 
since these clusters or supernebula are the sources of infrared, they are also the starburst. 

II Zw 40 is a starburst dwarf galaxy 10.5 (H Q /75) Mpc away. It appears to be the result 
of a collision of two smaller galaxies (Baldwin, Spinrad, & Terlevich 1982; Brinks & Klein 
1988; van Zee, Skillman, & Salzer 1998). At optical wavelengths it is dominated by one 
bright star cluster (Sargent 1970), which is also a strong source of Ha (Sargent & Searle 
1970) and associated nebular lines (Walsh & Roy 1993), has the Wolf-Rayet feature (Vacca 
& Conti 1992), and in general has all the signs of a young star formation region. The 
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infrared spectrum is dominated by high-excitation lines such as [Nelll] and [SIV] (Roche et 
al. 1991; Thornley et al. 2000; Beck et al. 2002). The radio continuum emission is nearly 
all thermal emission (Klein, Weiland, & Brinks 1991; Klein, Wielebinski, & Thuan 1984; 
Jaffe, Perola, & Tarenghi 1978; Deeg et al. 1993), and the emission is relatively compact 
(Wynn- Williams & Becklin 1986; Sramek & Weedman 1986; Klein, Weiland, & Brinks 
1991). Like the radio continuum, the Ha (Sargent & Searle 1970) and Brackett 7 lines 
(Wynn- Williams & Becklin 1986; Moorwood & Oliva 1988; Joy & Lester 1988; Vanzi et al. 
1996; Davies, Sugai, & Ward 1998) are bright in this galaxy . Molecular gas is conspicious 
by its absence; Sage et al. (1992) remark on its unusually high star formation efficiency 
In these respects, II Zw 40 resembles NGC 5253 (Crowther et al. 1999; Turner, Ho, & Beck 
1998; Turner, Beck, & Ho 2000), probably the youngest, most massive, and most extreme 
example of an obscured young super star cluster or supernebula in the local universe. 

Motivated by the galaxy's similarity to NGC 5253, we imaged II Zw 40 at 6, 3.6 
and 2 cm using the Very Large Array 7 in its A configuration for the maximum angular 
resolution, and at 11.7/im with the Long Wavelength Spectrometer (LWS) at the W. M. 
Keck Observatory. 8 The observations and results are described in the next section and the 
nature of the source in §3, in which we also compare the supernebula in II Zw 40 to other 
dwarf starbursts and discuss the significance of the structure seen in the radio maps. 

7 The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. 

8 The W. M. Keck Observatory is operated as a scientific partnership among the California 
Institute of Technology, the University of California, and the National Aeronautics and Space 
Administration. The Observatory was made possible by the generous financial support of 
the W. M. Keck Foundation. 
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2. Observations and Results 

2.1. Radio Images 

II Zw 40 was observed at 6, 3.6 and 2 cm on 23 August 1999, using the VLA in the A 
configuration. The observations began in late morning with reasonable summer weather. 
During the run the 12 GHz site monitor rms phase rose from 12 to 20°. Absolute fluxes 
were calibrated using 3C286, with fluxes of 3.45 Jy, 5.18 Jy, and 7.49 Jy at 2, 3.6, and 6 cm, 
respectively. We estimate absolute flux (scaling) uncertainties of ~ 5% at 6 and 3.6 cm, 
and ~7% at 2 cm due to source structure and variability in 3C286. Missing antennas in the 
inner array contributed to these uncertainties. For the 2 cm data, we used fast-switching 
with a 50s-calibrator/100s-source cycle time. We spent 6 minutes, realtime, on source 
at 6 cm and 3.6 cm, and 79 minutes, realtime, on source at 2 cm. However, given the 
fast-switching duty-cycle, the actual on-source time at 2 cm was about 44 minutes. The 
nearby calibrator 0552+032 was used for the fast- switching cycle; we measure a flux of 0.53 
Jy for this source. The fast- switching cycles were bracketed with standard calibration scans 
on the stronger calibrator 0532+075, for which we measure a flux of 1.65 Jy. We applied 
referenced pointing during the run, with pointing scans taken at 3.6 cm every hour on 
0532+075. In this way we were able to obtain reasonable images with the A configuration 
at 2 cm in daytime summer observing. We estimate, from our mapping of the calibrators, 
that the "seeing" at 2 cm with our fast-switching cycle is less than <0.01", which is the 
upper limit we fit to the deconvolved sizes of both calibrators. This is less than a tenth of a 
beam. Due to the lack of short spacings, the images are insensitive to structures that are 
more extended than ~ 30 times the beam. This approximate limit, 9 max , is listed for each 
wavelength in Table 1, along with beamsizes and rms noise levels. 

Images at each wavelength are in Figure 1. These maps were made with natural 
weighting. The emission is clearly resolved, especially in the highest resolution 2 cm map. 
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The position (J2000) of the 2 cm peak of 0.64 mJy/beam is a = 05 h 55 m 42:614 iCTOOl, 

5 = 03° 23' 32"03 ±0 / .'01. The main source is kidney-shaped, hinting that there is further 
unresolved substructure on sizescales less than CC/4 (20 pc). There is a secondary 2 cm 
peak, with a flux of 0.54 mJy/beam, located 0.23" (11 pc) to the east of the main source, 
at a = 05 h 55 m 42:629 ±0*001, 5 = 03° 23' 32"04 ±0"01. The 6 cm and 3.6 cm sources are 
consistent with the structure seen at 2 cm. 

Peak fluxes and other information for naturally- weighted maps are presented in Table 1. 
Total mapped fluxes were obtained by integrating the flux radially outward from the peak 
until it reached a maximum using AIPS task IRING. These mapped fluxes are presented 
in Figure 2 with the single dish fluxes and other interferometric fluxes from the literature. 
Most of the 6 cm and 3.6 cm flux is contained within a radius of 5-6". The 2 cm flux is 
confined to the inner 1"5 (75 pc) radius region. These sizescales are similar to max , the 
angular size at which undersampling effects set in. Sixty to eighty percent of the single dish 
emission is concentrated within this central few arcsecond portion of II Zw 40. 

The spectral index of radio continuum emission, which is defined as a in S u oc u a , is -0.1 
for pure thermal emission, some larger negative value, typically -0.6 to -1.2 for non-thermal 
emission from supernovae events and remnants, and takes on positive values when the 
radio emission is optically thick. We cannot formally derive the spectral index because 
the naturally-weighted beams at different wavelengths are not matched: the maps at the 
shorter wavelengths will lose relatively more flux due to undersampling. The spectral index 
of the overall radio emission in II Zw 40 has been discussed at length in Wynn- Williams 

6 Becklin (1986); Klein, Weiland, k Brinks (1991); Deeg et al. (1993). These studies 
indicate that the cm-wave emission of II Zw 40, even at 20 cm, is mostly thermal, and that 
thermal emission is overwhelmingly dominant at wavelengths less than 6 cm. 
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2.2. 12/xm Image 

II Zw 40 was observed on the night of 17 February 2000 using the LWS (Jones & 
Puetter 1993) on the Keck 1 telescope in imaging mode. The resulting image is shown in 
Figure Id (greyscale) and in Figure 3. The exposure was 486 seconds through a 1/rni wide 
filter centered at 11.7//m. (3 Gemini and a Bootes were used for calibration. The night was 
humid and non-photometric; scatter in the standards was 7% in an hour but as high as 
18% over the night, so we conservatively assign 15% uncertainty to the infrared flux. The 
pixel size is 0.08" and the chip 128 x 128; the actual resolution, based on the standard star 
images, is 0.3" to 0.5". The infrared source observed is either not or only slightly resolved; 
the FWHM is ~ 0.5". At this resolution the source cannot be distinguished from circularly 
symmetric. 

The total flux at 11.7/xm is 0.24 Jy. This agrees with the 0.22 Jy (Rieke & Low 1972) 
observed with a 6" aperture and is more than half of the 0.46 Jy seen by IRAS with a 
beam that took in the entire galaxy. The bright mid-infrared emission indicates that the 
compact radio source is not a supernova remnant. The II Zw \0 source has all the radio 
and infrared behavior of an HII region ionized by young OB stars. The ratio of mid-infrared 
to thermal 2 cm radio emission is ~ 20 for this compact radio/IR source, on the low side 
of the range expected for an HII region (Ho, Beck, & Turner 1990) but quite out of the 
range of supernovae. The value of 12yum/6 cm is only a factor of two greater than the value 
expected on the basis of Lya heating alone (Genzel et al. 1982), and it is nearly an order 
of magnitude smaller than it is in the otherwise very similar supernebula in NGC 5253 
(Gorjian, Turner, & Beck 2001). The presence of an infrared excess (IRE) over the value 
expected from Lya heating is generally attributed to factors such as the competition of 
dust and gas for UV photons and the absorption of longer wavelength light by dust. The 
conditions in II Zw 40 — low metallicity and a dominant OB population — may not be 
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conducive to the "normal" IRE. 

We list measured IR fluxes from the literature for II Zw 40 in Table 2, including 
large-aperture IRAS fluxes, and we plot the IR SED in Figure 2. We fit the SED of II Zw 40 
with three dust components, adopting dust emissivities proportional to v 13 , where (3 = 1.5. 
The fits are also illustrated in Figure 2. 

Our 11.7yum flux of 0.24 Jy agrees well with the fluxes measured by other groups (Rieke 
& Low 1972; Wynn- Williams & Becklin 1986) in few-arcsecond apertures, and suggests 
that the mid-IR flux of the central starburst region is confined to the 0.5" (25 pc) source 
in our image. The spectrum of this source is rather flat in the mid-IR. It therefore requires 
fairly warm emission; we can fit it with a dust component of temperature T~ 200K, giving 
a luminosity of L ~ 2.4 x 10 8 L for this component (see Figure 2). Dust emission from 
the supernebula in NGC 5253 and from the dwarf galaxy SBS 0335-052 also indicate hot 
(200K) dust, based on high resolution imaging at 11.7/xm and 18.7/xm (Gorjian, Turner, & 
Beck 2001; Dale et al. 2001). 

We can estimate a total Li R for II Zw 40 by including IRAS fluxes. The 25/iin IRAS 
flux of II Zw 40 is due to hot dust with a temperature of ~100K, we estimate, to explain 
both the 25^tm and 60/um components-again, like NGC 5253. This component (Figure 2) 
contributes about 7.2 x 10 8 L , about three times what the 200K dust contributes. The 
dust emission is unusually warm in this galaxy (Vader et al. 1993), although the cool 
(42K) dust does contribute about half (9.4 x 10 8 L ) of the total IR luminosity. The 100/xm 
is also complicated by the presence of Galactic cirrus. From the dust emission out to a 
wavelength of 100/xm we obtain a total infrared luminosity of L ~ 1.9 x 10 9 L . 
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2.3. Radio Analysis of the Source: Thermal and Nonthermal Emission and 

Sizescales 

There is a wealth of radio continuum observations of II Zw 40, both single dish and 
interferometric, and thus much information on the distribution of thermal and nonthermal 
emission in the source (Jaffe, Perola, & Tarenghi 1978; Klein, Wielebinski, & Thuan 1984; 
Sramek & Weedman 1986; Klein, Weiland, & Brinks 1991; Deeg et al. 1993).Vanzi et al. 
(1996) and Joy & Lester (1988) also discuss the radio structure and the thermal/nonthermal 
flux ratio. Since our A array images do not detect large spatial structures, it is useful 
to compare them to lower resolution maps and single dish fluxes which are sensitive to 
extended emission. 

At 6 cm, the lower resolution VLA images of Wynn- Williams & Becklin (1986) and 
Sramek & Weedman (1986) revealed 16 mJy of flux within a 20" region (corresponding 
to their undersampling radius). The total 6 cm single dish flux is about 21-22 mJy (Jaffe, 
Perola, & Tarenghi 1978; Klein, Weiland, & Brinks 1991). We detect 15±1 mJy at 
6 cm, 68% of the total single dish flux. Vanzi et al. (1996) predict, based on their 
extinction-corrected Ha flux, a 6 cm flux of 15.7 mJy for a 15" aperture. Joy & Lester 
(1988) also conclude that their near-IR and Br7 fluxes are in excellent agreement with the 
VLA fluxes. From this we conclude that at least 70% of the single-dish 6 cm emission is 
thermal, and confined to a region of less than 6" (300 pc) in extent. The remaining 30% of 
the single dish 6 cm emission forms a diffuse extended (>15", or 750 pc, based on the lower 
resolution VLA maps of Wynn- Williams & Becklin 1986) halo of thermal emission about 
this central, bright source. 

We detect a similar percentage of the total single dish emission at 2 cm. Our flux of 14 
mJy is ~75% of the estimated single dish emission of 18-19 mJy at 15 GHz (extrapolated 
from the 10.55, 10.7, and 24.5 GHz measurements of Jaffe, Perola, & Tarenghi 1978; Klein, 
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Weiland, & Brinks 1991; Klein, Wielebinski, & Thuan 1984). Our 2 cm flux is consistent 
with the VLA fluxes of Wynn- Williams & Becklin (1986) and Sramek & Weedman (1986) 
of 12 mJy from lower resolution observations. The 14 mJy of 2 cm flux is confined to within 
a 3" diameter region; of this, 10 mJy is found in the inner, 1.5" diameter, bright "core" 
region visible in the contour plot of Figure 1. The remaining 4-5 mJy of single dish flux 
that is not detected in any of the VLA maps must either be very diffuse, extended over 
sizescales >20" (1 kpc), or consist of a spatially extended collection of low flux sources. 

The total 2 cm flux is much greater than the sum of the two peaks seen in the 
naturally-weighted maps. For our flux of 10 mJy for the inner 1.5" (75 pc) core region, 
the required Lyman continuum fluxes are at least N Lyc = 1.0 x 10 53 s _1 , and higher if the 
emission has significant optical depth that it is partly thick at 2 cm. The total mapped 2 
cm flux of 14 mJy, which is contained within the inner 3", requires Nl yc = 1.4 x 10 53 s _1 . 

2.4. "Matched-Beam" Radio Images 

We cannot directly compare fluxes at different wavelengths for the naturally-weighted 
maps of the previous section. This is because the missing short baselines scale to different 
angular resolutions at the different wavelengths: the spatial scale at which emission is 
resolved out in the 2 cm maps (9 max ~ 3 — 4") is three times smaller than for the 6 cm map 
(#max ~ 10"). So the maps at each wavelength are missing a different fraction of the true 
total flux. 

It is possible, however, to make maps at different wavelengths for which fluxes can be 
directly compared. For this we must restrict the baselines in the inner (u,v) plane in the 6 
cm and 3.6 cm maps to match the 2 cm (u,v) restriction. With the current data set, we 
eliminated all baselines for which B/A < 20, 000. The resulting 6 cm and 3.6 cm images are 
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noiser than before, but their sensitivity to extended structure is close to that of the 2 cm 
image. We then convolved these "nearly identical" beams to the Gaussian beamsize of the 
6 cm map, 0"88 x 0"43, p. a. 52°, and compared them. 

Total mapped fluxes for these (u,v)-restricted maps, which resemble the naturally- 
weighted 6 cm map in appearance, are presented in the last column of Table 1. At 6 and 3.6 
cm the "matched beam" images contain less flux than do the naturally -weighted images: 
as the undersampling sizescale, 8 max , decreases, more of the flux is resolved out. Of the 15 
mJy of flux in the naturally- weighted 6 cm map (9 max ~ 10"), only 9 ±1.5 mJy remains in 
the "matched beam" image; of the 12 mJy of flux in the naturally-weighted 3.6 cm image, 
only 10±1.5 mJy is left in the "matched beam" map. But since the U band map is already 
restricted with respect to short spacings, the "matched beam" flux is the same as for the 
naturally-weighted map; the 2 cm flux is 14 mJy. 

Comparison of the "matched beam" fluxes indicates that the spectral index of the 
compact (9 < 4") emission in II Zw 40 is rising from 6 cm to 2 cm, with a value between 
+0.4 and +2. This strongly suggests that a significant fraction of the 2 cm free-free emission 
is not only compact, but optically thick. 

2.5. Highest Resolution Image: The Compactest Structure 

To further investigate the compact structure we have mapped the 2 cm emission with 
the highest possible resolution, and the result is presented in Figure 4. This map was made 
with uniform weighting, which places greater weight on the longer baselines. For standard 
VLA configurations, this sharpens the beam by roughly 50%. The noise in this image, 0.17 
mJy beam -1 , is a factor of two higher than in the naturally- weighted map, which is why 
less emission is seen. The uniform beam is 0.14"x 0.12", p. a. 44°. 
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The uniform image reveals 3-4 peaks within the central emission region, all of which 
are 0.59-0.60 mJy/beam, 3.5a. The three strongest peaks are separated by 05015, or 11 pc. 
They are at (J2000) a = 5 h 55 m 42 s .630, 42^615, 42*601 ± 0*001, 5=03° 23' 32'.'01, 32 / . , 02, 
32'.'02 ± O'.'Ol, from east to west, respectively. The first two sources listed also appear as the 
secondary and main peak in the lower- resolution 2 cm image of Figure 1. The fourth source 
is 6 pc to the north of the westernmost source of the triplet, also 6 mJy/bm, at a=5 h 55 m 
42!601 ± 0:001, 5=03° 23' 32'.'13 ± O'.'Ol. 

These sources may simply be the peaks of the more extended distribution, although 
there is an increase in brightness temperature with decreasing beam. The 2 cm naturally- 
weighted map has Tb ~ 120K, while the uniform map, with a beam a factor of 1.7 smaller 
in area, has T b ~ 200K. This suggests that there is indeed some unresolved structure. The 
brightest peak in the uniform map, at a = 5 /l 55 m 42*525, 5= 03° 23' 32'.'88 has a peak flux 
of 0.76 mJy/beam, ~ 5a. This peak is located about 1.5" to the northwest of the center of 
the triplet. There is little extended emission associated with this source, as it is not strong 
in the lower resolution images. 

We can estimate the source sizes using the spectrum and the intensities observed 
here. Peak flux densities of ~ 0.6 mJy/beam for the quadruplet correspond to brightness 
temperatures of ~200 K for the C'14 by / .'12 beam. In order for this brightness to be 
consistent with the rising spectrum, and optically thick emission, the source must be 
confined to a region smaller than the beamsize. For an electron temperature of 13000K, as 
determined from optical lines by Walsh & Roy (1993), and for r ~ 1, we estimate that a 
source diameter of d~ 1 pc would produce the observed brightness. The emission measure 
implied by a turnover frequency of 15 GHz is about 10 9 cm -6 pc, and the corresponding 
nebular density would then be ~ 3 — 4 x 10 4 cm -3 . Emission measures and densities of 
these magnitude are seen only in the youngest HII regions in the Galaxy. This suggests 
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that these nebulae are very young. The dynamical time for the expansion of an HII region 
to this size at a sound speed of 10 km/s is ~ 10 5 yrs; its size is another suggestion of its 
likely youth. 

For each source within the quadruplet, we find a lower limit to the Lyman continuum 
of N Lyc ~6x 10 51 s" 1 for t = (optically thin), N Lyc = 0.8 - 1 x 10 52 s" 1 for r = 1, and 
higher if r > 1. 

2.6. Comparison with Brackett 7; Extinction in the Nebula 

Brackett 7 fluxes have been obtained and analyzed by several groups, including 
Wynn-Williams & Becklin (1986) and Joy & Lester (1988). Vanzi et al. (1996) and 
Davies, Sugai, & Ward (1998) have actually imaged Brackett 7. We can compare our high 
resolution free-free fluxes with their results to estimate extinctions. 

We measure 14 mJy of free-free emission at 2 cm in II Zw 40, from a region ~3" 
in diameter. Using the recombination emissivities of (Hummer & Storey 1987) for 
12500K, we would predict on the basis of the 2 cm free-free flux a Br 7 flux of Sg^ = 
1.2 x 10~ 13 erg s~ x cm -2 . This can be compared to the Sg^ s 7 = 0.3 — 0.4 x 10~ 13 erg s~ x cm~ 2 
measured by Vanzi et al. (1996) in an aperture of 3-4", and S^ s 7 = 0.4 x 10~ 13 erg s _1 cm -2 
measured by Davies, Sugai, & Ward (1998). The extinction at 2.2 fim implied by the 
difference is A K ~ 1 — 1.2, for a visual extinction of A v ~ 8 — 10 mag (adopting the 
extinction curve of Rieke & Lebofsky 1985). Based on comparisons of Ha and Br7 
fluxes, Joy & Lester (1988), Vanzi et al. (1996), and Davies, Sugai, & Ward (1998) 
estimate A K ~ 0.32 — 0.36, a factor of 3 lower than we find. It is likely that the apparent 
contradiction is caused by the fact that much of the extinction is high and internal to the 
nebula. Infrared line ratios frequently give higher extinctions than do visible line ratios 
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(Beck et al. 1986). A K ~ 1 mag is seen in many starbursts of this luminosity (Ho, Beck, & 
Turner 1990; Kawara, Nishida, & Phillips 1989). 

We interpolate a single dish flux of 18.5 mJy at 2 cm from the 2.88 cm and 1.18 cm 
measurements of Klein, Weiland, & Brinks (1991) and Klein, Wielebinski, & Thuan (1984). 
For apertures larger than ~4" we would predict a Br7 flux of Sg r r c ^ = 1.6 x 1CT 13 erg s _1 cm~ 2 , 
which can be compared with the observed flux of = 7 x 1CT 14 erg s" 1 cm -2 measured 
by Vanzi et al. (1996) for a 15" aperture and Moorwood & Oliva (1988) in a 6" x 6" 
aperture. From these fluxes we estimate Ak ~ 0.96, or Ay ~ 8 for emission on larger 
(r > 200 pc) sizescales. 

3. The Nature and Stellar Population of the Source 

The infrared and radio emission from II Zw 40 has the form of a concentrated 1.5" 
(75 pc) diameter core containing at least half (10 mJy) of the large aperture flux, within 
a slightly larger (3", 150 pc diameter, 14 mJy) region of lower-level (2-2. 5a) extended 
emission, and diffuse thermal emission extended on larger scales which cannot be detected 
in these images. Our further discussion will concentrate on the compact thermal sources 
detected in these maps. 

The number of stars needed to ionize the radio source can be estimated from the total 
ionization of 1.4 x 10 53 sec -1 . This requires ~ 1.4 x 10 4 typical 07 stars over the central 
3" (150 pc) region, of which 10 4 are confined to the inner 75 pc core. If the stars follow the 
nebular emission — as seems to be the case in NGC 5253 — then from its appearance in the 
highest resolution map II Zw 40 is currently forming at least four pc-sized clusters with ~ 
600 O stars each. For a Salpeter ZAMS IMF down to 1 M Q , these four clusters would have 
masses of ~ 10 5 M & , and higher if the IMF extends to below 1 M . The estimated mass 
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for the entire 150 pc region would be at least ~2x 10 6 M (03 to G). 

The luminosity of this large population of young stars is consistent with the IR 
spectrum, in which the mid-infrared contribution is substantial. The total observed IR 
luminosity out to 100/rni is Li R ~1.9x 10 9 L . We estimate, using a ZAMS Salpeter IMF 
with upper mass cutoff of 03 (following Turner, Ho, & Beck 1998) that the ratio of ZAMS 
luminosity/Lyman continuum photons is NL yc ~ 2 x 10~ 44 s -1 . Then our Lyman continuum 
rate of 1.4 x 10 53 s _1 implies an OB luminosity of Lob ~ 3 x 10 9 L . This OB luminosity is 
very similar to Lob obtained by Vanzi et al. (1996) on the basis of emission line excitation. 
Lqb m &y exceed L IR , although there are uncertainties in the SED and in the fit, and in 
Lob- What is clear, however, is that the current starburst is responsible for the infrared 
luminosity; even a ZAMS OB luminosity appears to overpredict Lir. This suggests that the 
starburst is exceedingly young, or that there are unusually massive stars, 02 or above, as 
might be suggested by comparison with the smaller and less extreme cluster R136 (Massey 
& Hunter 1998). The thermal radio and infrared source, which we call the supernebula 
because of its similarity in all but size to a normal HII region, provides all of the infrared 
luminosity in II Zw 40. In fact, in this galaxy, the infrared luminosity may fall short of 
measuring the true star formation rate. 

We would like to emphasize a point made by Vanzi et al. (1996): the radio to IR 
correlation seen in many galaxies does not hold in II Zw 40, since the cm- wave fluxes fall 
short of the value predicted by the radio/IR correlation by factors of 6-10. Given that 
the IR excess in this galaxy appears to be low, based on the thermal radio fluxes, this 
discrepancy is even more notable. There is just not much nonthermal synchrotron emission 
in this BCD for the magnitude of its star formation. This is consistent with the starburst 
being extremely young. However, the very similar starburst in NGC 5253 shows a similar 
lack of nonthermal emission even though it has clusters estimated to range from 2.5-50 Myr 
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(Calzetti et al. 1997; Tremonti et al. 2001). It is unclear if the lack of nonthermal 
synchrotron emission from these starbursting dwarfs is due to the youth of the starburst or 
to some other cause, such as the lack of a global magnetic field (Turner, Ho, & Beck 1998). 

3.1. The Structure of the Starburst in II Zw 40: Multiple Star Clusters? 

The starburst in II Zw 40 clearly has a complex radio structure, best seen in our 2 cm 
map (Fig. la and Fig. 3). In the highest resolution uniform maps there are at least four 
sources visible, within a region of emission that extends to a diameter of 150 pc. Under 
the very rough assumptions of the earlier sections we estimate that the starburst contains 
14,000 07 star equivalents. The starburst appears to be made up of four compact, 1-pc 
diameter nebulae excited by clusters each containing ~600 O stars and ~ 10 5 M . The 
individual nebulae/clusters are separated by 10-12 pc. These characteristics are similar to 
the supernebulae in NGC 5253, He 2-10, and NGC 2146 (Beck et al. 1996; Turner, Ho, & 
Beck 1998; Kobulnicky k Johnson 1999; Turner, Beck, & Ho 2000; Tarchi et al. 2000; 
Beck, Turner, & Gorjian 2001). But it should be noted that what we refer to as "clusters" 
and "sources" are not isolated and well-defined emission regions, like globular clusters in 
our Galaxy. Rather, these are emission peaks within an extended region of emission. In 
this sense the nebulae resemble large Galactic HII region complexes such as W49 and W51. 
It would be very worthwhile to re-observe this galaxy with the highest possible resolution 
and better S/N. 

If the O stars deduced from the radio flux formed alongside the usual number of smaller 
stars, there will be a total mass of ~ 2 x 1O 6 M in young stars within the 150 pc volume 
of the starburst, and an estimated 10 6 — 10 7 young stars. This stellar density is high but 
not extraordinary; NGC 5253 for example has an estimated 1 million stars within a region 
1 by 2 pc (Turner, Beck, & Ho 2000). However II Zw 40 is a tiny galaxy. Sargent & Searle 
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(1970), Joy & Lester (1988) and Vanzi et al. (1996) all point out that the OB mass of 
II Zw 40 is a remarkably high (10%) fraction of its total stellar mass, which Vanzi et al. 
(1996) estimate at M ~ 2 x 10 7 M . 

The 75 pc extent of the starburst core places constraints on the trigger for the starburst. 
The crossing time for a disturbance moving at the sound speed in 10 4 K gas is roughly 
5-10 Myr. While it is difficult to place ages on these nebulae, comparison with Galactic HII 
regions of similar properties would give the nebulae ages of < 1 Myr (Gorjian, Turner, & 
Beck 2001). The ages of these sources are almost certainly a small fraction of the crossing 
time of the region, suggesting that the star formation was triggered nearly simultaneously. 
It seems unlikely, therefore, that "sequential super star cluster formation" is operating in 
II Zw 40. 

How does the structure seen in II Zw 40 compare to that in other supernebula- 
dominated systems? II Zw 40 is like NGC 5253 and He 2-10 in that the infrared and radio 
fluxes come from obscured super star clusters. The II Zw 40 source is more luminous than 
the NGC 5253 supernebula in the radio and infrared and contains at least twice the number 
of O stars, but they are distributed in multiple sources over a much larger volume. The 
total stellar content of the II Zw 40 source is close to that of the main source A in the 
dwarf galaxy He 2-10 (Kobulnicky & Johnson 1999; Beck, Turner, & Gorjian 2001; Vacca, 
Johnson, & Conti 2002), which it also resembles in size and radio and infrared flux. He 2-10 
has not been observed with high resolution so the possibility of discrete sub-clusters within 
the main sources is open. Why is the starburst in II Zw 40 concentrated within 75 pc, that 
in He 2-10 200 pc, that in SBS 0335-052 1000 pc (Dale et al. 2001; Izotov et al. 1997), 
and in NGC 5253 only 2 pc? That must depend on the interactions that triggered the 
formation of the massive star clusters in the first place and their subsequent evolution. 
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4. Conclusions 

Our high-resolution radio and infrared observations of the starburst dwarf galaxy 
II Zw 40 have found that the radio emission at 6, 3.6, and 2 cm is thermal free-free emission 
with ~70-75%, of the cm- wave emission confined to a region less than 3" in diameter. 
Within this region, there is a bright core of diameter ~ 1.5" (75 pc) containing at least 
four sub-sources. The spectrum of the compact emission appears to be rising between 3.6 
cm and 2 cm, suggesting that the free-free emission is optically thick at 2 cm, with implied 
emission measure of EM ~ 10 9 cm~ 6 pc and electron density n e ~ 3 — 4 x 10 4 cm -3 . For 
optically thick emission, the observed brightnesses require that the emitting sources are ~1 
pc diameter HII regions, with N Lyc ~6x 10 51 s _1 . The 75 pc "core" starburst region is a 
complex of four compact "supernebulae" , each powered by ~600 O stars, with an estimated 
~ 10 5 M Q in young stars in each cluster. The supernebulae/clusters are separated by 
~12 pc. The total ionization of II Zw 40 is at least N Lyc = 1.4 x 10 53 s _1 , or 14,000 07 star 
equivalents. For a Salpeter ZAMS IMF, we estimate that roughly 1-10 million young stars 
and a young stellar mass of 2 x 10 6 M are present within II Zw 40. The total young stellar 
luminosity is Lob ~3x 10 9 L Q . The OB luminosity is in good agreement with or perhaps 
even higher than the observed infrared luminosity of Lpr ~ 1.9 x 10 9 L . The radio-infrared 
ratios in this galaxy differ from the usual starburst values; we attribute this to the extreme 
youth of the star formation activity, the low metal content of the gas, and a very weak cool 
dust component. 
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Fig. 1. — Radio and infrared images of II Zw 40. Contours are ±2 n / 2 x 0.27 mJy/beam 
(~ 3er), starting at n=0. Beams are plotted in the lower right corners, a) top, left: 2 cm 
VLA map of II Zw 40. Beam is / /22 x / .'14, p. a. 46°. b) top, right 3.6 cm map. Beam is 
0"48 x 0"25, p.a. 52°. c) bottom, left: 6 cm map of II Zw 40. Beam is 0'.'89 x Of! A3, p.a. 52°. 
d) bottom, right: 6 cm contours atop the 11.7/im LWS image of II Zw 40. The seeing in the 
IR image is estimated to be 0^2 to C'5. 

Fig. 2. — Radio and infrared spectral energy distributions for II Zw 40. a) left: Radio SED 
of II Zw 40. Crosses with error bars: single dish measurements of Jaffe, Perola, & Tarenghi 
(1978), Klein, Wielebinski, & Thuan (1984), Klein, Weiland, & Brinks (1991), Deeg et al. 
(1993). Filled triangles: total fluxes from our naturally- weighted images. Open triangles: 
fluxes from the images with (u,v) range restricted to match the 2 cm (u,v) coverage (see text). 
Statistical uncertainties are ±1 mJy bm" 1 ; systematic uncertainties due to sidelobes are 
estimated to be 50% larger. The single dish data are fit with an optically thin bremsstrahlung 
spectrum, S oc v a , a — —0.1, with contributions of 21 mJy at 6 cm and 12 at 92 cm. The 
nonthermal synchrotron component is fit with a = —0.75. b) right: IR SED of II Zw 40. 
Fluxes and references are in Table 2. 

Fig. 3. — a) top, left: Optical image of II Zw 40 from Baldwin, Spinrad, & Terlevich (1982). 

b) top, right: The 11.7/im LWS image of II Zw 40. The entire frame is about 10" across. 

c) bottom, left: Blue optical image of II Zw 40 (Baldwin, Spinrad, & Terlevich 1982). d) 
bottom, right: LWS image placed atop the blue optical image to the same scale. The LWS 
inset is ~10" square. 

Fig. 4. — a) left: Uniformly-weighted 2 cm VLA map of II Zw 40, in greyscale, with naturally- 
weighted 2 cm map (see Figure 1) overlaid in contours. Greyscale range is 0.4-0.6 mJy, 
~2.5-4cr. Uniformly-weig hted beam is 0'.'14 x 0'.'12, p.a. 45°. Contours are 2 n/2 times ±0.4 
mJy/beam (~ 6a for the naturally-weighted map), starting at n=0. b) right: Uniformly- 
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weighted 2 cm map. Contours as for a), ~ 2.5a for this map. 



Table 1. II Zw 40 Radio Fluxes 



Single Dish r.m.s. Beam S h max Peak Flux Total Mapped Fraction of "Matched" Beam d 

Flux a Flux b - C Single Dish Flux Fluxes 

(mjy) (mjy/bm) "x", pa° (FWHM) (mjy/beam) (mjy) (mjy) (mjy) 



6 cm 22 0.11 0.89 x 0.43, 52 10" 3.4 15±1 0.68 9±1.5 

3.6 cm 21 0.088 0.48 x 0.25, 52 7" 1.4 12±1 0.57 10±1.5 

2 cm 18.5 0.10 0.21 x 0.14, 45 4" 0.64 14± 1.5 0.75 Mil. 5 

a Extrapolatcd from values at 4.9, 10.7, and 24.5 GHz. Klein, Wciland, & Brinks (1991), Klein, Wielebinski, & Thuan (1984), Jaflfc, Pcrola, & 
Tarenghi (1978) 

b 9max is the maximum sizescale that is well sampled by these images. Fluxes and peak fluxes are therefore lower limits to the total flux. 
c Total fluxes are obtained by integrating over a circular aperture centered on the peak flux in naturally-weighted maps. 

d "Matched" beam fluxes were obtained by making maps with (u,v) data restricted to baselines > 20kA; these maps are sensitive to structures 
smaller than 4" only. See text. 
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Table 2. Infrared Continuum Spectrum 





10.1 fim 


10.4 ^m 


11.7 fj,m 


12 jum 


20 fim 


25 jum 


60 /iin 


100 ^m 


Flux (Jy) 


0.20 a 


0.22 b 


0.24 c 


0.46 d 


1.0 e 


1.91 d 


QM d 


5.8 d 


Uncertainty (Jy) 


0.02 


0.04 


0.04 


0.05 


0.2 


0.2 


0.7 


0.9 



a Wynn- Williams & Becklin (1986) 
b Rieke & Low (1972) 
c This paper. 

d Vader et al. (1993). 100/xm flux is affected by Galactic cirrus emission. 
c Roche et al. (1991) 
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